Preprint typeset in JHEP style - HYPER VERSION 



Shear Dynamics in Bianchi I Cosmology 



Diego L. Caceres, Leonardo Castaneda and Juan M. Tejeiro 



Observatorio Astronomico Nacional, Universidad Nacional de Colombia, Bogota, 
Colombia 



E-mail: 


dlcaceresu@bt . unal . edu . co 


jmtejeirosOunal . edu. co 





Abstract: We present the exact equation for evolution of Bianchi I cosmological model, 
considering a non-tilted perfect fluid in a matter dominated universe. We use the definition 
of shear tensor and later we prove it is consistent with the evolution equation for shear 
tensor obtained from Ricci identities and widely known in literature ||, ||], |J. Our result 
is compared with the equation given by Ellis and van Elst in || and Tsagas, Challinor 
and Maartens Q. We consider that it is important to clarify the notation used in ||, || 
related with the covariant derivative and the behavior of the shear tensor. 
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1. 1+3 Orthonormal frame approach 

In a cosmological space-time (Ai, g) there are preferred worldlines representing the average 
motion of matter at each point, associated with comoving or fundamental observers, which 
do not have peculiar velocities. The signature used in this article is ( — h ++). The 4- 
velocity of the comoving particles is u a , u a = (1,0,0,0), u a u a = —1. This 4-velocity 
is orthogonal to the surfaces of spatial homogeneity. Therefore, it is defined the spatial 
projection tensor h a p as ||, [|l0|]: 



h a /3 := g a p + u a U/3, 



(1.1) 



given this tensor, we can define the orthogonally projected symmetric trace-free part of 
any tensor T a @ of second rank as: 



r (a/3> 
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rj-ryS 



X2) 



moreover, other two derivatives can be defined |3|: the covariant time derivative, along the 

75 ' 



fundamental worldlines, where for any tensor T a ^ • 



and the fully orthogonally projected covariant derivative V, where: 



(1.3) 



V e T^ = h<*X v h%h\h\V v T» v KX . (1.4) 

Given these derivatives, the 4-aceleration can be written as u a = u a .RU^. With these 
definitions the first covariant derivative of u a is decomposed into its irreducible parts, 
defined by their symmetry properties ||, [ 10 1 : 



V/3-U a = -U a U/3 + U a /3 + a a p + ^Qh aj3 , 



(1.5) 
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where a a p := u/ a .^)—Uf a u^) — ^Qh a p is the trace-free symmetric rate of shear tensor (o~ a p = 
°~{a0)i a a0U^ = 0, o~ a a = 0), which describes the rate of distortion of the matter flow; and 
uj a /s := U[ a .p] — U[ a up] is the skew-symmetric vorticity tensor (u a p = W[a/3], ^apu^ = 0), 
describing the rotation of the matter relative to a non-rotating (Fermi-propagated) frame 

§■ 

It is possible to obtain a propagation equation for the shear tensor, from Ricci identities 
§, §: 

_ U <V> = -le a -P + u^uP) - o-jfo-Ph - U {«J) - - \ a ^j , (1.6) 

where E a/3 is the Electric Weyl tensor, E a p = Ca^pgiflv? , where is the Weyl tensor, 

and 7T Q( g = Tj$h J ^ a h s ^, where T a p is the energy momentum tensor and ir a p is the trace-free 
anisotropic pressure. 

The Weyl Tensor is completely determined from its electric and magnetic parts, the last 
one defined as Q, Q defined as: 

H a p = ^VaSeC%u\ (1.7) 

where rj a ^ = u s r]g a p~ is a volume element for the rest spaces and r\ a p^s is the 4-dimensional 

volume element g {r] a p y s = V[a^8], %123 = y/\det(g a p)\). 

Using the Gauss-Codacci relation the Spatial Riemann Tensor is 0: 

3 R a /3 1 s = h> a h v ph e y h sR^e^ - VfU a Vsup + VgUa^-yU/3- (1.8) 
and the Spatial Ricci tensor is: 



2u - -9 2 + 2cr 2 + 2A 
P 3 



, (1.9) 



3 Ra/3 — ~^(a/3) ~ ® a af3 + V^'iig) + U^Up) + K a p + ~h a p 

where /i is the energy density and a 2 = \cF a ^a a p. 
2. Bianchi I cosmology 

Bianchi cosmologies are spatially homogeneous but not necessarily isotropic. For a review 
of Bianchi models, see [|J], 0], |J and for orthonormal frame approach Q, Q, ||. 

Here we will consider Bianchi I cosmology. The metric of this model is given by Q, ||, 

ds 2 = -dt 2 + a 2 (t)dx 2 + b 2 {t)dy 2 + c 2 (t)dz 2 , (2.1) 
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and the average expansion scale factor S(t) = y/ abc. It reduces to the FLRW case when 
a(t) = b(t) = c(t) = S(t). Given this metric the connection components are: 



r* x — ad, r^ t — - , 



r* = bb T\ = - (2-2) 

1 yy UU, L yt , V J 



T* = rr Y z 
x zz 00 j x zt 



We are going to study the dynamic evolution of shear tensor from these connection com- 
ponents. 

2.1 Shear Dynamics 

The solution for the scale factors can be obtained directly from Einstein Equations when 
we consider a perfect fluid Q: 

a(i) = 5(t)exp(EiW(t)), (2.3) 
b(t) = S(t)expp 2 )W(t), (2.4) 
c(t) = S(t)eM^sW(t)), (2.5) 

where 

W(t) = f-JL, (2-6) 
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and the constants T, a satisfy Q: 
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From (|J), (U) and we obtain: 

T .... = V, 

-S»(t) 



(2.10) 



o^x = El > (2-11) 
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a. 



c 2 (t) 



(2.12) 



(2.13) 



J3 'S 3 (t)' 

Now, the spatial homogeneity of the Bianchi I space-times ensures that all invariants depend 
at most on time. It is an irrotational universe, u a p = 0, and also it is spatially flat, 
S Ra/3 = 0. With a perfect fluid, = 0, the equation Ql.Sj ) can be reduced to pi, 0: 



cr a p 



(2.14) 



Following this equation, it seems that Van Elst and Ellis in Q and Tsagas, Challinor and 
Maartens in || conclude that in the absence of anisotropic pressures the shear behaves as 
S- 3 . From ( gig) it is easy to give a non-correct interpretation for the shear dynamics 
because we could conclude shear tensor is equal to a constant times S~ 3 , but from ( 2.11| ), 
( |2~l2|) and (|2~lg| ) we see the scale factors play a role in shear dynamics. 
For checking our result we verify our shear expression (pTTlp , (|2l2| ) and ( |2~13D is consistent 



with shear evolution equation (1.6). Given the definition of covariant derivative we have: 



da-. 



dt 



a 
a 



(2.15) 



It is not immediate to integrate this equation to get a a S~ s (t). Taking into account this 
term in the evolution equation: 

da TT (a b c\ . 

(2.16) 



On 



dt \a b c I 
Now, we will see our shear satisfies the evolution equation for the shear tensor: 



da,, 



dt 




(2.17) 



We check it satisfies the evolution equation (1.6). 



Now, we present the generalized Friedmann equation. It is an equation that allows us to 
integrate the scale factor S(t). We assume a 7-law for state equation (p = (7 — 



M 



(2.18) 
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where [i is the energy density and M = 0, the dot the covariant time derivative defined in 
fll.3j ). So, it can be shown || that the generalized Friedmann equation is given by: 



s 2 _ r 2 m 

where 2T 2 = V 2 + r 2 , + r 2 . When T = we get the usual Friedmann equation. The 
additional term at right is the shear Energy. 



For late times the shear constant T does not play a significant role in the evolution of scale 
factors, but at early times it has a great difference with the FLRW model. It can play 
an important role in physics processes in the early universe, such as Nucleosynthesis and 



structure formation ||, il|| , [12] , [jjjf| . For dust, we get an analytic solution for S(t): 



S(t) = y ^Mt 2 + y/SEt. (2.20) 

It is a different expression that the one shown in Q, it can be checked it satisfies the 
Friedmann equation ( |2.19| ). Using this expression for S(t) we get: 



= — ln( , , 4S , (2.21) 



\t + 



s/3M , 



Thus, with these analytic solutions it is straightforward to obtain a(t),b(t) and c(t), and 
therefore the components of the shear tensor. We consider M = 1 and three cases for the 
constants Si, E2 and S3: 



1. Si = 0, S 2 = -S 3 = 0.05. 

2. Si = S 2 = -0.05, S 3 = 0.1. 

3. Si = S 2 = 0.05, S 3 = -0.1. 

When it is considered the limit t — >• 0, S(t) — >■ 0, there are two types of singularities, 
the cigar singularity, which is the case 1 and 3 and the pancake singularity which is the 
second case @], 0. The cigar case means that two of the scale factors tend to while the 
third increases withouth bound, while the pancake case means that one of the scale factors 
tend to and the other two increase. 



In figures 1, 2 and 3 we present the evolution of Shear components for these three cases. 



We conclude these components tend to zero, as we could expect from the equation for shear 
tensor we got before. However, the appropiate way to know if this model isotropize is to 
consider the evolution of two scalars S 2 and W 2 P] defined as: 
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Figure 1: Shear Evolution for case Si = 0, £2 = — £3 
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Figure 2: Shear Evolution for case Si = £2 



0.045 



-0.05, £3 = 0.1. 



2 = (2>22) 

W = — ( 2 - 23 ) 

These quantities are defined because the components of shear tensor are not dimension- 
less, they are normalized with the Hubble scalar H, and hence measuring the dynamical 
importance of the different variables with respect to the overall expansion of the universe. 

When both parameters tend to zero we can say the model tends to isotropy. As was pointed 
out by H , it was thought Bianchi VIIq non-tilted dust model isotropize in the sense S — >• 
when t — > 00. However, for t — > 00 W — > Wo, where Wo is a constant whose value can 
be any positive number depending on the initial conditions [||. On the contrary these two 
factors tend to zero for our three cases of Bianchi I model. We illustrate these behaviours 
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Figure 3: Shear Evolution for case Si = £2 = 0.05, £3 
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Figure 4: Anisotropy parameters E 2 and W 2 for case Si = 0, £2 



0.05 



in figures 4, 5 and 6. 



We see these parameters tend to zero in the three cases we have considered, the model 
isotropize. In the first case, as the constant Si = 0, the shear component a xx = 0. 
However, the electric component E xx , which we have not plotted here but we showed in 
1 14], is not identically zero, although it is very small, compared with the other Electric 
components. Both Electric and Shear tensors are diagonal. In the other cases, where there 
exists axial simmetry, it is reflected in the components of shear tensor. 



3. Conclusions 

We have shown the shear tensor in BI cosmology and we analized the solutions in the 
dust model. Prom our analysis it is clear is convenient to be careful when the covariant 
derivative is considered. 




The anisotropy parameters show us that this model isotropize, for different cases. For late 
times S(t)a t 2 / 3 , Ha t~\ a afi a afi a S(t)~ 6 a t~ 4 , so Y?a t~ 2 , while E^E^a S(t)~ 12 , so 
W 2 a t~ 8 and as we see in the plots, these parameters decay, given a well defined behavior 
of the kinematical quantities in BI cosmology. 
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